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Abstract 


Experiments  were  performed  on  two  simulated  porous 
media  made  up  of  different  packing  arrangements  of  0.5  inch 
lucite  spheres  between  two  lucite  plates.  A  Newtonian  oil 
was  used  as  the  fluid  medium. 

A  laser  doppler  velocimeter  was  chosen  to  monitor  the 
velocity  in  the  pore  space  and  to  determine  the  onset  of 
turbulence.  Velocity  measurements  in  the  pore  space  were 
made  by  detecting  the  doppler  frequency  shift  of  a  laser 
light  which  was  scattered  by  contaminant  particles  in  the 
oi  1 . 

The  study  successfully  demonstrated  the  use  of  the 
velocimeter  for  interstitial  velocity  measurements  over  a 
range  of  flowrates.  The  results  indicated  that  a  porous 
medium  should  not  be  represented  by  a  series  of  straight 
flow  channels  of  varying  cross  sections.  The  flow  geometry 
should  include  a  curvilinear  path  of  the  fluid. 
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1.  Introduction 


Pioneering  work  in  the  study  of  flow  through  porous 
media  ,  particularly  in  the  turbulent  region,  was  initiated 
at  the  University  of  Alberta  by  McFarland  (1).  His  objective 
was  to  see  if  the  onset  of  turbulence  in  a  porous  medium 
occured  at  a  single  value  of  Reynolds  number.  The  method  of 
streaming  bi ref r i ngence  was  used  for  flow  visualisation.  The 
experiment  was  limited  in  that  the  fluid  used  was 
non-Newtonian,  and  there  was  insufficient  pressure  head  to 
creat  turbulence.  This  work  is  a  continuation  of  that  effort 
using  the  method  of  laser  doppler  velocimetry.  Like  the 
method  of  birefringence,  this  method  is  also  non- i ntrusi ve . 
However,  a  Newtonian  fluid  can  be  used  with  this  technique. 

The  Darcy  equation  has  long  been  used  for  predicting 
the  flow  of  fluids  through  porous  media.  It  has  been 
recognized  that  above  a  certain  flowrate,  the  pressure  drop 
across  the  medium  no  longer  increases  linearly  with  the 
flowrate,  though  the  flow  may  still  be  laminar.  The 
deviation  from  the  linear  relationship  is  due  to  the 
increasing  importance  of  inertial  effects  in  the  porous 
medium.  Ultimately  the  flow  becomes  turbulent  as  the  flow 
velocity  increases. 

Throughout  the  literature  several  flow  equations  have 
been  suggested  to  describe  non-Darcy  flow.  Of  these,  the 
Forchheimer  equation  has  been  the  most  well  received.  It  has 
been  shown  to  have  a  mechanistic  basis  by  Irmay  (2).  The 
Reynolds  number  has  been  commonly  used  to  define  the  limits 
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of  validity  of  the  flow  equations.  This  study  proposes  to 
use  a  Newtonian  fluid  and  to  define  the  upper  limit  of  the 
laminar  region  by  a  Reynolds  number.  The  second  objective  is 
to  investigate  if  turbulent  flow  beyond  could  be  described 
mathematical ly . 

The.  measuring  technique  adopted  for  this  study  is  laser 
doppler  velocimetry  operating  in  the  fringe  mode.  The  method 
is  unique  in  that  it  is  possible  to  obtain  local  fluid 
velocity  measurement  in  undisturbed  porous  media.  The 
velocimeter  records  the  instantaneous  velocity  of  tiny 
impurities  flowing  with  the  fluid.  When  these  naturally 
occur ing  particles  are  small  enough,  it  is  reasonable  to 
assume  that  they  travel  at  the  same  velocity  as  the  fluid.  A 
profile  of  the  local  velocity  provides  information 
concerning  the  nature  of  the  flow  in  the  interstitial 
regions  of  the  porous  medium.  Streamline  or  laminar  flow 
would  be  indicated  by  a  steady  velocity  profile  over  a 
period  of  time.  The  onset  of  turbulence  is  character ized  by 
a  continuous  random  fluctuation  of  velocity.  This  condition 
would  then  give  rise  to  an  unsteady  velocity  profile  which 
changes  with  time. 
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2.  Review  of  Literature 


Numerous  efforts  have  been  made  to  understand  the  flow 
of  fluids  through  porous  media,  particularly  in  the  region 
where  the  Darcy  equation  is  no  longer  applicable.  Early 
observers  attributed  the  non-Darcy  flow  to  the  onset  of 
turbulence,  assuming  a  similarity  to  the  flow  in  a  straight 
pipe.  Hubbert  (3)  argued  that  comparing  flow  in  a  porous 
medium  to  that  in  a  straight  pipe  is  not  justified  as  in  a 
porous  medium,  non-Darcy  flow  appears  well  before  turbulence 
sets  in.  This  has  been  verified  experimentally  by  several 
researchers  (4, 5, 6, 7, 8). 


2 . 1  Flow  Equations 

In  the  case  of  a  one  dimensional,  horizontal  Darcy 
flow,  the  pressure  gradient  across  the  medium  is  directly 
proportional  to  the  flowrate  through  it,  i.e. 


(1) 


where  ^  is  the  pressure  gradient,  u  is  viscosity  of  the 
fluid,  K  is  the  permeability  of  the  the  medium  and  q  is  the 
superficial  velocity.  To  describe  the  flow  beyond  the  Darcy 
region,  Forchheimer  (9)  in  1901  presented  the  following  two 
equations : 


3 


4 


(2) 


and 
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where  a,  b,  c  and  m  are  constants.  He  suggested  the  use  of 
Equation  2  as  i t  could  be  reduced  to  the  Darcy  equation  at 
low  values  of  q.  Since  then  Equation  2  has  been  called  the 
Forchheimer  equation.  Other  researchers  (10)  continued  to 
use  Equation  3  and  values  of  m  between  1  and  2  have  been 
suggested.  This  equation  is  inadequate  in  that  the  constants 
do  not  in  any  way  describe  the  physical  characteristics  of 
the  porous  medium.  Though  it  is  easier  to  fit  Equation  3  on 
log- log  paper,  Trol lpe  et  al.  (11)  have  noted  that  with  more 
reliable  experimental  data,  a  Forchheimer  relation  can  be 
fitted  more  accurately.  In  general,  Equation  2  is  more 
widely  accepted.  It  can  be  derived  from  the  Navier-Stokes 
Equation  as  shown  by  Irmay  (12)  and  Chywl  (13).  For  a  fully 
saturated  homogeneous  linear  horizontal  system,  Chywl  has 
shown  that  Equation  2  takes  the  form  : 
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where  p  is  the  density  of  the  fluid,  and  Fa  and  Fb  are 
constants  characterizing  the  porous  medium.  They  are  often 
called  the  viscous  and  inertial  resistance  coefficients, 
respectively.  Comparison  of  Equation  4  and  the  Darcy 
equation  suggests  that  Fa  is  the  inverse  of  the  permeability 

a 

of  the  medium.  Thus  Equation  4  may  be  written  as 

‘  S7  =  I  q  +  Fbpq2  (5) 


2.2  Definition  of  Reynolds  Number 

The  Reynolds  number  is  often  used  to  define  the  limits 
of  the  validity  of  the  various  flow  equations.  A  common  form 
of  the  Reynolds  number  used  in  the  literature  has  been 


N  =  Dpa  (6) 

r  y 

where  D  is  the  average  particle  or  sphere  diameter.  In  some 
cases,  D  has  been  taken  to  be  the  average  diameter  of  the 
pore  space.  The  disadvantage  of  using  Equation  6  is  that 
neither  of  the  definition  of  D  may  be  a  character istic 
length  of  the  porous  medium.  Green  and  Duwez  (14)  presented 
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another  expression  which  could  be  obtained  from  the 
definition  of  the  Reynolds  number,  which  is  the  ratio  of  the 
inertial  to  the  viscous  forces.  Comparing  the  correspondi ng 
terms  in  Equation  5,  the  Reynolds  number  may  be  written  as 


This  expression  is  preferred  as  it  has  a  mechanistic  base. 
It  has  been  derived  directly  from  the  definition  of  the 
Reynolds  number.  The  expression  has  the  length  parameter 
KF^  ,  which  is  a  product  of  the  two  characteristics  of  the 
porous  medium.  Thus  KF^  may  be  a  better  characteristic 
length  than  D. 


2.3  Correlation  Parameter 

A  typical  correlation  used  in  flow  studies  through 
pipes  and  porous  material  is  the  relationship  between 
friction  factor  and  Reynolds  number.  Using  this  correlation, 
Chalmers  et  al.  obtained  results  similar  to  those  obtained 
for  coiled  tubes  (15).  This  similarity  suggested  to  them  the 
possibility  of  approximating  flow  through  porous  media  to 
the  flow  through  coiled  pipes.  Fancher  and  Lewis  (16) 
plotted  Fanning  friction  factor  versus  and  obtained  a 
separate  curve  for  each  porous  medium.  They  suggested  that 
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such  plots  could  be  used  to  indicate  the  gradual  change  from 
viscous  to  turbulent  flow  taking  place. 

Green  and  Duwez  (17)  defined  the  friction  factor  as  the 
ratio  of  'dissipative'  forces  to  inertial  forces.  Using  this 
definition  and  Equation  5,  the  friction  factor  is  given  as 


AP 


(8) 


With  the  definition  of  N  from  Equation  7,  Equation  8  may 

1  r 

then  be  writ  ten  as 


C 


f 


+  1 


(9) 


When  Cf  was  plotted  against  N  ,  Green  and  Duwez  (18) 

r  p 

showed  that  all  the  data  from  the  various  porous  media  fell 

along  the  curve  given  by  Equation  9.  As  the  choice  of  N 

for  the  Reynolds  number  is  preferred,  plots  of  Cf  versus 

N  wi 1 1  be  used  in  this  study, 
rp 
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2.4  Visualization  Techniques 

The  dye  tracer  technique  is  a  popular  method  used  in 
many  early  flow  studies.  Among  those  who  used  it  were 
Schneebeli  (19),  Chauvateau  and  Thirriot  (20),  Wegner, 
Karabelas  and  Hanratty  (21)  and  Kyle  and  Perrine  (22). 

Though  the  method  is  simple,  its  main  disadvantage  is  that 
the  coloured  streamlines  tend  to  break  up  after  moving  a 
short  distance  downstream.  Furthermore,  it  will  contaminate 
a  closed  system. 

Streaming  bi ref r i ngence  is  another  proven  technique  in 
flow  visualisation.  It  is  an  excellent  qualitative  method 
for  the  study  of  laminar,  transition  and  turbulent  flows. 
McFarland  (23)  used  this  method  in  his  studies. 

Unfortunately  pure  liquids  exhibit  bi ref r i ngence  only  at 
very  high  shear  rates.  Flow  at  moderate  rates  would  then 
require  the  use  of  colloidal  suspensions  which  are 
non- Newton i an . 

The  hot-wire  anemometer  gained  popularity  in  dynamic 
fluid  studies  in  the  1960's.  Wright  (24)  was  one  of  the 
first  to  adapt  the  technique  for  use  in  a  porous  medium. 
Since  interstitial  velocity  can  be  measured,  the  anemometer 
is  useful  for  detecting  any  form  of  velocity  fluctuation  in 
the  pore  region.  The  onset  of  turbulence  could  then  be 
noted.  Kingston  and  Nunge  (25)  and  van  der  Merwe  and  Gauvin 
(26)  also  used  the  hot-wire  anemometer  in  their  experiments. 

The  main  disadvantage  of  hot-wire  anemometry  is  that 
the  technique  is  intrusive.  Its  very  presence  alters  the 
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conditions  in  the  flow  stream.  The  hot-wire  can  be  made  very 
small.  It  is  very  fragile  and  has  to  be  carefully  inserted 
in  the  pore  space. 

Laser  doppler  velocimetry  appears  to  be  an  attractive 
technique  for  quantitative  flow  studies  especially  in  the 
turbulent  region.  The  phenomenon  of  doppler  frequency  shift 
has  been  Known  for  many  years.  It  was  not  until  1964  when 
Yeh  and  Cummins  (27)  successfully  demonstrated  the 
application  of  laser  doppler  velocimetry.  They  conducted 
measurement  of  a  steady  laminar  flow  of  water  in  a  pipe 
using  a  He-Ne  laser.  Their  success  prompted  others  to  use 
the  technique  (28,29,30).  Goldstein  and  Kreid  (31)  reported 
accuracy  of  measurement  of  0.1%.  Measurements  in  the 
turbulent  region  were  first  made  by  Goldstein  and  Hagan  (32) 
in  1967. 

This  technique  proved  popular  and  by  the  early  1970's, 
complete  commercial  packages  were  available  for  velocity 
measurement  (33).  The  main  advantage  of  this  method  is  that 
it  is  non- i n t rus i ve  and  thus  does  not  perturb  the  flow.  It 
also  records  velocity  directly  and  requires  no  calibration 
(34).  However,  in  this  method  the  medium  must  be  optically 
transparent  to  the  laser  beam.  There  must  be  sufficient 
scattering  particles  in  the  fluid  to  cause  a  strong  doppler 
signal.  In  some  cases,  these  scattering  particles  may  have 
to  be  introduced  into  the  fluid. 

To  date,  only  one  study  of  flow  in  a  porous  medium 
using  a  laser  doppler  velocimeter  has  been  reported.  In 
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1974,  Johnston  et.  al.  (35)  measured  velocities  in  the 
interstitial  regions  of  a  packed  cell  made  up  of  transparent 
Plexiglas  spheres.  The  fluid  used  had  a  refractive  index 
close  to  that  of  the  glass  spheres  at  the  wavelength  of  the 
laser  beam.  The  investigation  demonstrated  that  a  laser 
velocimeter  could  be  used  for  flow  measurements  in  a  porous 
med i urn . 


3.  Theory  of  Measurement 


3.1  Principle 

A  laser  doppler  velocimeter  is  an  optical  flowmeter. 

The  light  beam  from  a  laser  source  is  scattered  by  particles 
moving  with  the  flowing  fluid.  A  doppler  frequency  shift 
occurs  and  it  is  defined  as 


f 


f 

s 


(10) 


where  f$  and  f.j  are  the  frequencies  of  the  scattered  and 
incident  beams,  respectively.  Written  in  terms  of  the 
velocity  of  the  moving  particle,  Equation  10  has  been  shown 
by  Watrasiewicz  and  Rudd  (36)  to  be 


(11) 


where  U  is  the  velocity  of  the  particle,  A  is  the  wavelength 
of  the  incident  beam,  e$  and  e^  are  the  unit  vectors  of  the 
scattered  and  incident  beams,  respectively.  From  Figure  1, 
it  can  be  deduced  that  the  velocity  component  measured  is  in 
the  direction  (  e  -  e .  ) .  The  magnitude  of  (  e  -  e-  )  can 
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Fig  .  1  DOPPLER  SHIFT 


be  shown  to  be  2. sin  (9/2),  where  9  is  the  angle  through 
which  the  incident  beam  has  been  scattered.  Rewriting 
Equation  11  in  terms  of  the  velocity  gives 


13 
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(12) 


Equation  12  shows  that  the  doppler  shift  f^  is  directly 
proportional  to  the  velocity  of  the  particle.  The  small 
particle  is  assumed  to  have  the  same  velocity  as  the  flowing 
fluid. 

The  variation  of  velocity  of  flow,  given  by  Equation 
12,  would  with  time  indicate  the  Kind  of  flow  behaviour 
taking  place  in  the  region  under  study.  A  constant  velocity 
is  characteristic  of  streamline  or  laminar  flow.  A  velocity 
profile  randomly  fluctuating  with  time  would  indicate 
turbulent  flow.  The  transition  between  the  laminar  and 
turbulent  regimes  has  sometimes  been  referred  to  as  a  period 
of  intermittent  flow.  In  this  region,  the  flow  alternates 
between  the  laminar  and  the  turbulent  states  (37). 


3.2  Mode  of  Operation 

Several  optical  arrangements  are  available  for  laser 
doppler  velocimeters .  They  are  the  ;  (1)  Reference  Beam  mode 
(2)  Differential  Doppler  mode  and  (3)  Dual  Scatter  Mode.  A 


■ 
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description  of  these  modes  of  operation  is  given  by  Durst  et 
a  1 . (38)  and  Watrasiewicz  and  Rudd  (39).  They  are  generally 
described  in  terms  of  the  doppler  effect.  However,  in  1969, 
Rudd  (40)  proposed  an  alternative  model  which  gives  a  good 
visual  description  of  the  system.  He  described  the  effect  in 
terms  of  a  fringe  pattern,  thus  the  name,  fringe  model. 


3.3  Fringe  Model 

A  beam  splitter  placed  in  front  of  the  laser  unit 
produces  two  coherent  monochromatic  laser  beams.  When  these 
two  beams  interfere,  they  form  a  set  of  dark  and  bright 
fringes  in  the  intersecting  volume.  The  spacing  between  the 
alternate  bright  fringes  (41)  is  given  as 


A  particle  moving  across  the  interference  fringes  with  a 
velocity  component  perpendicular  to  the  fringe  planes  will 
vary  the  light  intensity  of  the  fringe  pattern.  When 
observed  from  any  direction,  this  variation  is  in  the  form 
of  a  sinusoidal  intensity  fluctuation.  A  photomultiplier  is 
used  to  detect  this  fluctuation.  The  output  of  the 
photomultiplier  is  an  electrical  signal  having  a  frequency 
dependent  upon  the  velocity  of  the  particle  moving  across 


. 
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the  fringe  pattern  (  See  Figure  2  ).  The  frequency  of  light 
variation  may  be  written  as 


f 


m 


_1 

At 


(  14) 


where  At  is  the  time  to  cross  one  fringe  spacing.  It  may  be 
written  as 


At 


As  _  X/(2. sine/2) 
U  U 


(15) 


Therefore 
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U  .  2  .  s  i  n  9  /  2 

X 


(  16) 


Equation  16  is  identical  to  Equation  12  obtained  from  the 
doppler  shift  consideration. 

For  measurement  in  a  liquid,  Equation  16  becomes 


U  .  2  .  s  i  n 


V2 


(  17) 
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Fig.  2  FRINGE  SPACING 
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where  9^  is  the  angle  between  the  intersecting  beams  and  A, 
is  the  wavelength  in  the  liquid.  From  the  definition  of 
refractive  index  and  the  laws  of  refraction,  Longhurst  (42) 
showed  that 


A 


l 


n 


( 18) 


and 


n^sin(0£/2)  =  n  sin(e/2) 


(  19) 


where  is  the  refractive  index  of  the  fluid.  Therefore 


sinfQ^/z) 


—  s  i  n  ( 0  /  2  ) 


(20) 


or 


si n (0£/2 ) 


s  i  n  (  0  /  2  ) 

A 


(21  ) 


18 


Substituting  Equation  21  into  Equation  17  gives  Equation  16. 
Thus  Equation  16  is  also  valid  for  measurement  in  liquid 
provided  the  proper  values  of  6  and  X  are  used. 

3.4  Measuring  Volume 

The  volume  formed  by  the  two  intersecting  beams, 
sometimes  called  the  measuring  or  control  volume,  is  an 
important  parameter  in  the  velocimetry  system.  It  defines 
the  spatial  resolution  of  the  system. 

The  beam  intersecting  angle,  9,  is  calculated  from  the 
geometry  formed  by  the  two  beams  leaving  the  beam  splitter 
(  See  Appendix  1  ).  From  a  consideration  of  the  Gaussian 
properties  of  the  light  beams  and  the  optics  of  the 
converging  lens,  the  dimensions  of  the  measuring  volume,  as 
indicated  in  Figure  3,  has  been  shown  (43)  to  be 


(22) 


.  4  .  (X)(f)  .  _ 1 


D  cosQ/2 

e 


(23) 


Fig.  3  DIMENSIONS  OF  MEASURING  VOLUME 
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where  f  is  the  focal  length  of  the  converging  lens  and  De 
is  the  beam  diameter.  Using  Equations  13  and  23,  the  number 
of  fringes  in  the  measuring  volume  may  be  estimated  as 


(24) 


3.5  Turbulent  Flow 

The  instantaneous  velocity  at  any  point  in  the  flow  may 
be  written  in  the  form 


U  =  Urn  +  u 


(25) 


where  Urn  is  the  time  averaged  mean  velocity  and  u  is  the 
fluctuation  about  the  mean.  For  steady  laminar  flow,  there 
is  no  fluctuation  of  velocity  and  so  u  =  0  (44)  and  U  and  Urn 
are  equal.  In  turbulent  flow,  there  is  a  fluctuation  of  the 
instantaneous  velocity  about  the  mean  Urn.  However,  the  time 
averaged  value  of  this  fluctuation,  u  is  0.  As  mentioned 
earlier,  the  intermittent  region  is  character ized  by 
alternating  periods  of  laminar  and  turbulent  flow.  The 
fraction  of  time  during  which  turbulence  occurs  at  a  given 
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position  is  defined  by  Schlichting  (45) 
factor,  y.  Thus  y  =  0  defines  the  flow 
corresponds  to  the  onset  of  turbulence, 
intermittent  flow,  y  has  a  value  greater 
than  1 . 


as  the  i ntermi t tency 
as  1  ami nar  and  y  =  1 
During  the  period  of 
than  0  and  less 


4.  Experimental  Equipment 


4.1  Porous  Medium 

To  simulate  porous  media,  packed  cells  of  different 
flow  geometry  were  used.  It  was  originally  intended  that  all 
the  4  cells  from  McFarland's  work  (46)  were  to  be  tested  in 
this  study.  However,  only  two  of  the  cells  could  be  located. 
After  some  preliminary  runs,  it  was  found  that  only  one  of 
them  had  pore  spaces  large  enough  to  accommodate  the  two 
intersecting  laser  beams  without  reflection  from  the  sides 
of  the  spheres.  This  was  the  Cubic  No.  1  cell  with  spheres 
of  0.50  ±  0.01  inch  diameter.  The  stainless  steel  ball 
bearings  were  replaced  by  lucite  spheres.  Another  cell  of  a 
different  pack  was  constructed.  It  will  be  referred  to  as 
the  Cubic  No.  2  cell.  It  has  the  same  size  spheres  as  Cubic 
No.  1.  The  cubic  unit  cell  is  similar  to  that  in  Cubic  No.  1 
except  that  it  has  been  rotated  through  45s  about  one  of  its 
edges.  The  two  packing  arrangements  are  illustrated  in 
Figure  4  and  Figure  5. 

Each  pack  consisted  of  two  layers  of  spheres  between  1 
inch  thick  lucite  plates.  The  spheres  in  each  layer  were 
recessed  to  half  their  diameters  into  each  lucite  plate  to 
reduce  wall  effects.  Each  recess  in  the  lucite  was  made  with 
the  aid  of  a  programmable  drill  press,  producing  a  grid  of 
holes  of  centres  0.50  ±  0.005  inch  apart.  Both  cells  had 
arrays  of  10  rows  by  27  columns  of  spheres.  Each  cell  had  an 


22 


■ 

' 


23 


CUBIC 


No.  1 


CUBIC  No.  2 


Fig. 4  UNIT  CELLS  OF  PACKING  ARRANGEMENT 
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CUBIC  No.  1 


CUBIC  No.  2 


Fig.  5 


VIEW  OF  PARTIAL  CROSS  SECTION  OF  CELLS 
NORMAL  TO  FLOW 
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inlet  flow  diffuser  to  minimize  flow  perturbation.  Pressure 
taps  were  located  eight  sphere  columns  downstream  of  the 
inlet  and  four  sphere  columns  upstream  of  the  outlet.  The 
physical  properties  of  each  cell  are  presented  in  Table  1. 


4.2  Flow  System 

In  order  to  obtain  a  constant  pressure  head  sufficient 
to  cause  turbulent  flow,  oil  was  gravity-fed  from  a  tank  on 
the  eighth  floor  of  the  Chemi ca 1 -Mi nera 1  Engineering 
building.  A  schematic  of  the  flow  system  is  illustrated  in 
Figure  6.  A  continual  overflow  was  provided  at  the  upper 
tank  to  maintain  the  constant  head  of  pressure.  A  45-meter 
head  was  available  at  the  basement  of  the  same  building.  It 
was  a  closed  system  with  a  centrifugal  pump  to  return  the 
oil  to  the  upper  tank.  The  pump  used  was  an  Aurora  Turbine 
pump  driven  by  a  single-phase  1.5  hp  motor.  The  plumbing 
consisted  mainly  of  PVC  components.  The  two  tanks  were  made 
of  polyethylene.  The  properties  of  the  oil  used  in  the 
system  are  shown  in  Table  2. 

An  inlet  manifold  was  designed  such  that  three 
flowmeters  were  connected  in  parallel,  each  having  its  own 
control  valve.  Two  of  the  meters  were  Rotameters  while  the 
third  was  a  Trident  volumetric  flowmeter.  Depending  on  the 
flowrate,  the  two  Rotameters  could  be  used  in  parallel. 

The  two  pressure  taps  of  the  cells  were  connected  to  a 
Va 1 i dyne  differential  pressure  transducer.  Two  diaphragms  of 
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Table  1.  Physical  Properties  of  Packed  Cells 


Cell  # 
No . 


1 


(  cm  ) 

A  (  cm 2  ) 

Ax  (  cm 

1  .27 

16.42 

20.38 

1  .27 

11.88 

12.60 

2 
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TANK 


Fig.  6  SCHEMATIC  OF  FLOW  SYSTEM 


Table  2.  Properties  of  the  Oil  at  29.0C  ±  1.0‘C. 


V i scos i ty 


16.4 


mPa .  s 


Dens i ty 


850 . 1  Kg/m 3 
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maximum  rating  of  6.89  kPa  (  1.0  psi  )  and  0.69  KPa  (  0.1 
psi  )  respectively  were  available.  The  sensitivity  of  the 
system  was  1  %  of  the  maximum  rating  of  the  diaphragm.  This 
arrangement  measured  the  pressure  drop  across  the  two  taps. 


4,3  Laser  Equipment 

The  velocimetry  system  employed  in  this  study  used  the 

fringe  mode  with  forward  scatter.  The  laser  and  its 

associated  equipment  are  illustrated  in  Figure  7.  The  He-Ne 

laser  unit,  Spectra  Physics  (  Model  120  ),  produces  5  mW  of 

0 

continuous  luminous  radiation  of  a  wavelength  of  6328  A.  The 
Optical  Unit  split  the  incoming  beam  into  two  beams  which 
intersected  in  the  region  of  interest.  The  light  scattered 
by  the  moving  fluid  was  detected  by  the  photomultiplier.  An 
oscilloscope  was  necessary  to  monitor  the  presence  of  a 
doppler  signal.  For  this  purpose  the  Hewlett  Packard  564  was 
used.  The  doppler  frequency  signal  was  then  processed  and 
displayed.  The  display  could  either  be  in  the  form  of  a 
direct  frequency  read-out  from  a  meter  or  a  trace  on  an  X-Y 
recorder.  In  the  latter,  the  root -mean-square  value  of  the 
voltage  represented  the  fluctuation  of  the  measured  velocity 
about  the  mean.  The  velocimeter  used  in  this  study  is  the 
Type  55L  Laser  Doppler  Anemometer  from  Disa  Electronics.  The 
components  of  the  package  are  listed  in  Table  3. 
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LASER  VELOCIMETRY  SYSTEM 


Table  3.  Components  of  Type  55L  Laser  Doppler  Anemometer 


Laser  Unit  Spectral  Physics  Model  120 

Optical  Unit  Type  55L01 

Photomultiplier  Type  55L10 

PM  Tube  EMI  Type  9658B 

HV  Supply  Type  5 5 L 1 5 

Signal  Processor  Type  55L20 


32 


4.4  Mounting  of  Cell 

As  Goldstein  and  Kreid  (47)  had  observed,  it  was  more 
convenient  to  have  the  test  cell  mounted  on  a  moveable 
platform.  This  allowed  different  pore  sections  to  be 
examined  without  upsetting  the  laser  and  the  photomultiplier 
tube  alignment.  For  this  purpose,  a  milling  machine  was  used 
on  which  the  cell  rested  on  a  steel  plate  placed  on  three 
pneumatic  dampers.  The  milling  table  was  restricted  to 
movement  in  the  plane  normal  to  the  direction  of  the  laser 
beams.  Both  the  laser  and  optical  units  rode  on  an  optical 
bench  on  one  side  of  the  cell.  The  photomultiplier  and  the 
processing  electronics  were  on  the  other  side  of  the  cell. 
Plate  1  shows  a  view  of  the  optical  arrangement  and  the 
mounting  of  the  cell. 
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PLATE  1.  OPTICAL  ARRANGEMENT  AND  CELL  MOUNTING 


5.  Procedure 


5.1  Flowing  Medium 

Initial  tests  were  carried  out  using  tap  water  as  the 
flowing  medium.  The  pressure  drop  across  the  measuring 
points  was  too  low  to  be  recorded  by  the  differential 
pressure  transducer.  The  flowing  medium  was  then  changed  to 
a  blend  of  two  oils.  Uniflo  Base  Oil  and  Light  Atmospheric 
Gas  Oil  were  mixed  in  a  proportion  of  about  3  to  2  by 
volume.  The  Newtonian  nature  of  the  mixture  was  verified  on 
a  Weissenberg  Rheogoniometer .  It  was  unnecessary  to  seed  the 
oil  as  there  were  sufficient  scatterers  in  the  system  to 
provide  good  doppler  signals. 

It  was  observed  during  preliminary  runs  that  the  oil 
flowing  through  the  system  reached  a  steady  temperature  of 
29.0  (  ±  1.0  )*C.  All  subsequent  observation  and  measurement 
of  properties  were  made  at  that  temperature. 

The  properties  of  the  oil  are  given  in  Table  2.  The 
viscosity  of  the  oil  was  measured  with  the  aid  of  a 
Cannon- FensKe  viscometer  in  a  bath  Kept  at  29.0  C.  The  oil 
density  was  determined  using  a  density  bottle. 
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5.2  Calibration 

The  Rotameters  were  calibrated  by  collecting  volumes  of 
oil  passing  through  in  a  given  period  of  time.  The  amount  of 
oil  collected  was  weighed.  Calibration  graphs  of  mass 
flowrate  versus  Rotameter  reading  were  then  obtained.  These 
are  presented  in  Figures  A.1  and  A. 2. 

The  Trident  volumetric  flow  meter  was  checked  for  its 
accuracy.  A  volume  of  0.20  cu.  ft.  of  oil  as  recorded  on  the 
meter  was  collected  and  weighed.  The  volume  of  the  amount  of 
oil  collected  was  then  calculated  from  its  density  and 
weight.  The  mean  value  was  found  to  be  5670  ±  10  cc.  The 
flowrate  could  then  be  calculated  by  recording  the  time 
taken  for  0.20  cu.  ft.  (  meter  reading  )  of  oil  to  pass 
through  the  meter. 

The  Validyne  differential  pressure  transducers  were 
calibrated  in  the  Standards  Laboratory  of  the  Department  of 
Chemical  Engineering.  The  associated  equipment  consisted  of 
a  demodulator  and  a  digital  voltmeter.  Transducer  2  could 
not  be  calibrated  for  pressure  differences  less  than  0.19 
kPa.  The  regulator  valve  used  in  the  Standards  Laboratory 
could  not  be  finely  controlled  at  those  low  values.  The 
calibration  graphs  of  the  two  transducers  are  found  in 
Figures  A. 3  and  A. 4. 
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5.3  A1 iqnment 

Prior  to  any  optical  alignment,  the  tables  on  which  the 
laser  unit  and  the  packed  cell  were  to  be  placed  were 
isolated  from  the  floor  to  minimize  building  vibration. 

The  beam  from  the  laser  was  made  to  pass  through  a  beam 
splitter.  The  two  beams  were  then  aligned  to  ensure  that 
they  intersected  at  the  focal  point  of  the  convex  lens  of  13 
cm  focal  length.  Fine  adjustments  were  made  with  the  aid  of 
a  projecting  lens  positioned  at  the  beam  intersection.  With 
this  arrangement,  interference  fringes  could  be  projected  on 
to  a  screen.  The  beam  splitter  was  adjusted  until  the  best 
set  of  fringes  was  obtained.  After  the  proper  adjustments, 
the  projecting  lens  was  removed.  The  cell  was  symmetrically 
positioned  such  that  the  intersecting  volume  was  at  the 
center  of  the  flow  channel  in  a  pore  space.  Its  external 
surfaces  were  thoroughly  cleaned  to  ensure  minimal  optical 
distortion  of  the  laser  beams.  The  measuring  point  was 
arbitrarily  chosen  in  the  region  between  the  two  pressure 
taps.  Subsequent  runs  were  conducted  at  the  same  point.  The 
photomultiplier,  with  a  convex  collecting  lens  of  20  cm 
focal  length,  was  set  on  an  optical  bench.  For  optimal 
signal  strength,  it  was  positioned  symmetrically  between  the 
two  diverging  beams.  Using  the  inspection  window  at  the  side 
of  the  photomultiplier  assembly,  a  focussed  image  of  the 
beam  intersection  was  obtained.  The  pin-hole  aperture  of 
diameter  0.1  mm  was  adjusted  to  fall  on  the  intersecting 
volume.  This  step  was  repeated  several  times  until  the  best 
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s i gna 1  - to-noi se  ratio  was  obtained.  Plate  2  shows  a 
photograph  of  the  beam  intersection  (  a  red  X  )  located  in  a 
pore  space  of  cell  2 . 

With  a  doppler  signal  on  the  oscilloscope,  the  high 
voltage  supply  to  the  photomultiplier  was  set  to  provide  an 
anode  current  of  5  micro-amps.  The  amplifier  section  of  the 
tracker  was  adjusted  to  obtain  the  best  signal  waveform  on 
the  oscilloscope.  The  doppler  frequency  could  be  read  off 
directly  from  an  analogue  meter.  The  value  was  checked 
against  the  frequency  estimated  from  the  oscilloscope 
waveform.  This  was  to  ensure  that  the  tracker  had  locked  on 
to  the  correct  doppler  signal.  Another  output  from  the 
frequency  tracker  was  connected  to  the  X-Y  recorder  to  trace 
the  velocity  profile.  It  was  driven  by  its  internal  time 
base  of  2  sec/inch.  The  detailed  adjustment  procedure  of  the 
tracker  electronics  can  be  found  in  the  manufacturer's 
handbook  for  the  velocimeter  (48). 

A  needle  valve  at  the  inlet  to  the  cell  was  originally 
intended  for  controlling  the  flow  through  the  system.  It  was 
observed  that  controlling  the  flow  from  this  point 
introduced  flow  perturbations  too  close  to  the  pressure 
measuring  points.  The  gate  valves  of  the  Rotameters  were 
found  to  provide  better  control  of  the  flow. 

At  each  flowrate,  the  pressure  drop  across  the  taps  was 
recorded.  The  time  averaged  local  velocity  was  recorded  on 
the  X-Y  plotter.  Runs  were  repeated  to  ensure 
reproducibi 1 i ty .  The  flow  to  cell  2  was  reversed.  In  this 
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PLATE  2. 


VIEW  OF  INTERSECTING  BEAMS 
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case,  the  cell  was  referred  to  as  cell  2A.  This  was 
undertaken  to  see  if  the  location  of  the  pressure  taps  was 
critical.  The  pressure  transducer  of  6.89  kPa  rating  was 
used  for  cell  1.  Lower  pressure  differences  were  encountered 
in  cells  2  and  2A  and  the  transducer  of  0.689  kPa  rating  was 
used . 

It  must  be  mentioned  that  the  initial  alignment  was 
time  consuming.  However,  once  completed,  it  was  necessary 
only  to  check  for  the  focussed  image  of  the  intersecting 
volume  seen  through  the  inspection  window  of  the 
photomultiplier,  and  that  the  pin-hole  aperture  was 
positioned  over  the  image  before  each  run. 

The  beam  intersecting  angle,  0,  used  in  this  study  was 
the  largest  possible  with  the  Optical  Unit.  A  large  angle 
was  necessary  to  reduce  the  size  of  the  measuring  volume. 

For  a  more  accurate  study  of  the  local  velocity,  the 
measuring  volume  must  be  as  small  as  possible.  This  would 
require  a  larger  intersecting  angle  to  be  obtained  using 
some  other  optical  arrangement.  The  size  of  the  measuring 
volume  in  relation  to  the  cross  section  of  the  pore  space  is 
presented  in  Appendix  1. 


5.4  Treatment  of  Data 

After  the  data  had  been  processed,  a  linear  least 
squares  fit  of  the  Forchheimer  equation  was  carried  out  on 
each  set  of  data  to  obtain  estimates  of  k  and  F.  .  As  the 
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validity  of  the  Forchheimer  equation  in  the  turbulent  region 
has  not  been  established,  only  data  points  up  to  the  onset 
of  turbulence  were  considered. 

Inspite  of  the  damping  effects  of  the  recording  pen, 
the  X-Y  recorder  was  used  to  chart  the  velocity  profiles.  It 
was  chosen  as  it  was  convenient  for  recording  a  continuous 
profile  over  a  period  of  time  and  it  presented  sufficient 
qualitative  results.  For  more  detailed  and  quantitative 
results,  the  output  of  the  frequency  tracker  should  be  fed 
to  some  other  electronic  recording  devices  such  as  a 
magnetic  tape  recorder . 
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6.  Results 


The  optical  characteristics  of  the  system  are  presented 
in  Table  4.  The  calculation  of  these  properties  is  found  in 
Appendix  1.  A  large  intersection  angle  was  necessary  to 
reduce  the  size  of  the  measuring  volume.  However,  the  size 
of  the  angle  was  limited  by  the  maximum  setting  of  the  beam 
splitter  and  the  clearance  between  two  adjacent  spheres  in 
the  cell.  It  can  be  observed  from  Equation  16  that  a  larger 
angle  would  give  rise  to  a  higher  doppler  frequency. 

No  measurements  of  the  porosities  of  the  cells  were 
undertaken.  The  values  indicated  in  Table  6  were  calculated 
from  the  physical  dimensions  of  each  cell. 

The  experimental  and  the  processed  data  are  tabulated 
in  Appendix  3  in  Tables  A.1  through  A. 6.  A  sample 
calculation  may  be  found  in  Appendix  2.  The  graphs  showing 
the  relationship  between  the  pressure  drop,  A  P,  and  the 
flowrate,  Q,  are  found  in  Figures  8  through  10. 

Figures  11  through  13  contain  the  velocity  profiles 
recorded  on  the  X~Y  plotter.  The  vertical  axis  of  each  of 
the  traces  represents  the  doppler  frequency  or  indirectly 
the  instantaneous  local  flow  velocity  in  the  pore  space.  The 
horizontal  axes  represent  the  time  base  of  2  sec/inch.  All 
the  traces  have  been  presented  on  a  scale  of  1  :  1. 

A  steady  trace  indicates  a  steady  local  velocity  or 
streamline  flow  in  the  laminar  region.  At  this  point,  the 
i n termi t tency  factor,  y,  is  0.  As  the  flowrate  increases, 
per i od i c  ve 1 oc i ty  f 1 uctua t i on  i s  observed  (  See  Figure  11  ). 
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Table  4.  Properties  of  Optical  Intersection 


Beam  intersecting  angle  21.4* 

Width  of  intersecting  volume  0.869  mm 

Diameter  of  intersecting  volume  0.164  mm 

Fringe  spacing  0.0017  mm 

No.  of  interference  fringes  96 
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This  is  then  interpreted  as  the  start  of  the  intermediate 
zone  between  the  laminar  and  the  turbulent  states. 

Turbulence  is  reached  when  the  velocity  becomes  completely 
unsteady  in  time  and  y  is  equal  to  1.  For  an  accurate 
determination  of  the  onset  of  turbulence,  a  numerical 
analysis  of  the  instantaneous  velocity  is  necessary. 

However,  a  reasonable  estimate  of  the  onset  of  turbulence  is 
summarized  in  Table  5. 

Plate  3  shows  the  oscilloscope  waveforms  of  some 
typical  bursts  of  doppler  signals.  The  upper  picture  shows  a 
burst  of  doppler  signal  caused  by  the  presence  of  a  particle 
in  the  measuring  volume.  The  lower  picture  shows  another 
burst  of  signal  as  a  particle  was  leaving  and  another  was 
entering  the  measuring  volume.  The  arrival  of  particles  in 
the  measuring  volume  was  completely  random.  Noise  was  always 
present  in  the  system  and  it  could  be  differentiated  from 
the  doppler  signal  by  its  lower  amplitude  and  irregular 
periodicity.  A  reasonable  signal  strength  of  0.5  volt 
peak-to-peak  was  obtained. 

The  pressure  gradient  versus  average  velocity 
relationships  are  presented  in  Figures  14  through  16.  These 
graphs  also  contain  the  least  squares  fits  of  the 
Forchheimer  equation.  Table  6  shows  the  values  of  k  and  Fb 
obtained  for  each  cell  from  the  least  squares  fit. 

The  plot  of  the  friction  factor,  Cf,  versus  Reynolds 
number,  Nrp  is  displayed  in  Figure  17.  The  plot  contains 
data  from  Cells  1  and  2.  It  shows  the  curve  as  predicted  by 


« 
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Equation  9.  The  onset  of  turbulence  for  each  cell  has  also 
been  indicated.  Figure  18  is  a  similar  plot  containing  only 
data  from  Cells  2  and  2A. 
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PRESSURE  DROP  VS  VOLUME  FLOWRATE  DATA  FOR  CELL  t*2 
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FIG.  10  PRESSURE  DROP  VS  VOLUME  FLOWRATE  DATA  FOR  CELL  »2A 
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Q  =  50.4  cc/s  f,  =  139  kHz  N  =  0.035 

d  rp 


Q  =  92.5  cc/s  f,  =  268  kHz  N  =  0.064 
'  d  rp 


50  kHz/in. 


Q  =  186.9  cc/s  fd  =  495  kHz  Nrp  =  0.129 


Q  =  245.2  cc/s  fd  =  680  kHz  =  0.170 

150  kHz/in.  2  sec/in. 


FIG.  11 


VELOCITY  WAVEFORMS  FOR  CELL  1 
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Q  =  258.6  cc/s  f,  =  720  kHz  N  =  0.179 

a  rp 


Q  =  284.6  cc/s  f,  =  735  kHz  N  =  0.197 

'  d  rp 


Q  =  298.1  cc/s  =  760  kHz  N  =  0.206 


i  i 


150  kHz/in. 


2  sec/in. 

- *■ 


FIG.  11  -  CONTD. 
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Q  =  16.47  cc/s  f,  =  33  kHz  N  =  0.135 

'  d  rp 


Q  =  32.35  cc/s 


fJ  =  48  kHz  N  =  0.266 
d  rp 


Q  =  65.29  cc/s 


=  49  kHz 
d 


N  =0.536 
rp 


Q  =  67.29  cc/s  fd  =  47  kHz  Nrp  =  0.553 


15  kHz/in. 


2  sec/in. 


FIG.  12 


VELOCITY  WAVEFORMS  FOR  CELL  2 
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FIG.  12  -  CONTD. 
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Q  =  19.53  cc/s  fd  =  36.5  kHz  N  =  0.163 


Q  =  44.58  cc/s  f  =  57  kHz  N  =  0.372 

a  ro 


Q  =  61.05  cc/s 


fd  =  54  kHz  N  =  0.510 


Q  =  67.05  cc/s  f,  =  53.5  kHz  N  =  0.560 

'  d  rp 

« i 


15  kHz/in. 


2  sec/in. 


FIG.  13 


VELOCITY  WAVEFORMS  FOR  CELL  2A 
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Q  =  76.23  cc/s  f,  =  52.5  kHz  N  =  0.637 

d  rp 


Q  =  33.0  cc/s  f,  =  56  kHz  N  =  0.777 

d  rp 


15  kHz/in. 


2  sec/in. 


FIG.  13  -  CONTD. 
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Table  5 .  Critical  Values  at  Onset  of  Turbulence 


Cell  Q  f  N 

d  r 

#  (cc/s)  KHz 

1  258.6  720  104 

2  75.29  45.5  41.7 

2A  76.23  52.5  42.2 


N 

r  p 

0.  179 

0.618 

0.627 
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Y  axis: 

0.2  vo It /cm 

X  axis 

:  0.2  ms/cm 

Cel  1  1 

:  .  Q  =  18.0 

cc/s  f 

=42.5  KHz 

PLATE  3. 


OSCILLOSCOPE  WAVEFORMS  OF  DOPPLER  SIGNALS 
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PRESSURE  GRAD ] ENT  VS  SUPERFICIAL  VELOCITY  -  CELL  #1 
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Table  6.  Character i st ics  of  Porous  Media 


Cell  §  Porosity  K  ( x 1 0“ 7  ) 

%  m2  1  /  m 

1  47.64  4.7395  48.60 

2  47.64  4.8313  389.5 

2A  47.64  4.9967  383.0 
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REYNOLDS  NO  .  ,  Nrp 

FIG.  18  PLOT  OF  FRICTION  FACTOR  VS  REYNOLDS  NO.  -  CELL 


7.  Discussion 


A  small  scatter  of  data  is  observed  in  both  cells  with 
less  appearing  in  Cell  2.  In  the  case  of  Cell  2,  values  of 
AP  were  recorded.  This  was  due  to  the  shorter  distance 
between  the  pressure  taps  in  Cell  2.  Pressure  differences 
less  than  0.19  kPa  have  been  extrapolated  from  the 
calibration  of  pressure  transducer  2. 

The  data  of  Cell  1  were  recorded  with  the  intersecting 
volume  positioned  symmetrically  in  the  pore  space.  Runs  were 
repeated  with  the  intersecting  volume  placed  as  close  as 
possible  to  the  smallest  regions  of  the  pore  space.  Random 
velocity  fluctuations  were  observed  at  lower  flowrates  than 
with  the  intersecting  volume  at  the  widest  section  of  the 
pore  space.  However,  McFarland  (49)  reported  uniform 
occurrence  of  turbulence  in  the  pore  space.  The 
bi ref r i ngence  technique  is  a  qualitative  method  and  its 
visual  quality  depends  on  the  magnification  of  the  optical 
system.  It  is  possible  that  the  velocity  fluctuation  in  the 
smaller  sections  of  the  pore  space  could  not  be  observed 
with  the  visual  resolution  he  had  in  his  study. 

Flow  through  the  narrower  sections  of  the  pore  is 
faster  and  can  give  rise  to  waKes  and  vortices.  The  vel oc i t y 
fluctuations  in  these  sections  recorded  by  the  velocimeter 
may  have  been  due  to  these  f low  phenomena. 

In  Cell  2,  turbulence  occured  almost  at  the  same 
superficial  velocity  in  most  parLS  of  the  pore  space.  The 
data  reported  were  those  with  the  measuring  volume  placed  in 
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the  middle  of  the  pore  space. 

An  inspection  of  the  cross  section  of  Cell  1  (  see 
Figure  5  )  which  has  a  converging  and  diverging  cross 
section  suggests  that  the  flow  would  be  mainly  rectilinear. 
On  the  other  hand,  curvilinear  flow  would  be  more  likely  to 
occur  in  Cell  2.  This  was  verified  when  streams  of  air 
bubbles  were  occasionally  observed  to  travel  downstream  in 
the  cells.  This  observation  is  interesting  and  does  help  to 
explain  the  following  results. 

Figure  14  shows  that  with  Cell  1,  the  Forchneimer 
equation  fitted  with  data  points  up  to  the  onset  of 
turbulence  could  not  be  used  to  predict  the  pressure  drop  in 
the  turbulent  region.  On  the  other  hand  with  Cell  2,  Figures 
15  and  16  indicate  that  the  Forchheimer  equation  provides 
reasonable  prediction  of  the  pressure  drop  even  for  flow 
well  into  turbulence. 

Another  interesting  observation  is  illustrated  in 
Figure  17  which  shows  a  plot  of  the  friction  factor,  Cf, 
versus  Reynolds  number,  Nfp.  It  has  been  discussed  earlier 
in  a  previous  section  that  the  use  of  the  N^p  is  preferred 
over  N  for  the  definition  of  the  Reynolds  number.  As  such, 
only  plots  of  C^r  versus  N^p  will  be  shown.  Such  plots 
demonstrate  the  agreement  between  the  experimental  data  and 

the  theory  defining  C^  and  Nrp. 

After  the  onset  of  turbulence  in  Cell  1,  there  is  a 
departure  of  the  data  points  from  the  theoretical  curve 
predicted  by  Equation  9.  This  resembles  a  simi iar 
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correlation  in  the  case  of  a  straight  pipe.  There  is  good 
agreement  between  the  data  and  the  theory  in  the  case  of 
Cell  2.  Unfortunately,  a  comparison  of  the  same  plot  for 
Cell  1  by  McFarland  (50)  cannot  be  made  as  he  could  not 
record  any  observation  of  turbulence  with  the  glycerine 
solution  due  to  the  small  pressure  head  available  (51). 

The  results  shown  in  Figure  17  suggest  that  it  may  be 
inaccurate  to  characterize  a  porous  medium  as  a  series  of 
straight  flow  channels  with  converging  and  diverging  cross 
sections.  Houpert  (52)  considered  a  porous  medium  to  have  a 
a  cross  section  similar  to  that  of  cell  1  in  his  derivation 
of  the  quadratic  form  of  the  pressure  gradient  and 
superficial  velocity  relationship.  Though  his  equation  may 
have  the  form  of  the  Forchheimer  equation,  his  approximation 
of  the  porous  medium  is  questionable  since  the  flow  would  be 
mainly  rectilinear.  Houpert  (53)  and  Katz  et  al .  (54) 
attribute  the  extra  pressure  loss  in  non-Darcy  flow  to  the 
inertial  effects  of  the  fluid  flowing  through  constricting 
passages.  Though  this  may  be  partly  the  case,  it  is 
reasonable  to  assume  that  curvilinear  flow  makes  a  greater 
contribution  to  the  extra  pressure  drop.  Jones  (55) 
interprets  the  departure  from  Darcy  flow  as  being  caused  by 
secondary  flow  due  to  the  curvature  of  the  pores.  Such  flows 
have  been  observed  in  early  experiments  with  curved  pipes  by 

White  (56)  and  Taylor  (57). 

The  velocity  traces,  Figures  11  through  13,  indicate 
that  the  intermediate  region  in  Cell  1  extends  over  a  wider 
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range  of  flowraie.  This  intermediate  region  is  characterized 
by  increasing  periods  of  'spiky'  velocity  fluctuations. 

These  observations  are  similar  to  those  referred  to  by 
Schlichting  (58)  in  his  discussion  of  the  transition  zone  in 
a  straight  pipe.  No  clearly  defined  intermediate  zone  was 
observed  with  Cell  2.  A  higher  magnitude  of  the  fluctuation, 
u,  was  observed  in  Cell  1  than  in  Cell  2.  This  indicates 
that  there  was  higher  intensity  of  turbulence  in  Cell  1. 

The  effect  of  flow  geometry  is  also  evident  in  the 
results  of  Table  4.  In  particular,  Cell  2  has  a  much  higher 
inertial  resistance  coefficient  Fb  ,  although  both  cells  have 
the  same  value  of  porosity.  This  suggests  that  Fb  is 
dependent  on  path  curvature.  It  can  be  observed  from  the 
processed  data  in  Appendix  3  that  for  a  particular 
superficial  flow  rate,  a  higher  pressure  gradient  occurred 
in  Cell  2  than  in  Cell  1.  The  higher  values  of  Fb  and 
pressure  gradient  suggest  that  the  flow  in  Cell  2  was  more 
curvilinear  than  in  Cell  1.  Table  7  compares  the  results  of 
Cell  1  with  those  obtained  by  McFarland  (59).  There  is  a  one 
and  a  half  times  difference  between  the  estimates  of  Fb 
obtained  from  the  two  studies.  The  repack  of  the  cell  could 
have  contributed  to  the  difference.  The  more  probable  cause 
may  be  the  fact  that  there  were  insufficient  data  points  in 
McFarland's  estimation  of  Fb  from  the  data  using  the 
glycerine  solution. 

Table  8  contains  a  comparison  of  the  Reynolds  number,  Nr 
at  the  onset  of  turbulence  with  those  of  some  authors  in  the 


' 


•  ; 


66 


literature.  The  values  of  Nr  from  this  work  are  104  and  41.9 
(an  average  of  the  values  for  Cells  2  and  2A  ).  It  is 
observed  that  the  results  of  Cell  2  show  better  agreement 
with  those  of  Schneebeli  (60)  and  Kyle  and  Perrine  (61). 

Both  groups  of  researchers  conducted  experiments  in  which 
the  flows  were  mainly  curvilinear. 

The  results  obtained  here  support  the  view  that  it  is 
reasonable  to  use  kF^  as  the  length  parameter  in  the 
definition  of  Reynolds  number  although  no  single  value  of 
Reynolds  number  could  be  found  to  describe  the  onset  of 
turbulence.  The  use  of  the  Cf  vs  Nrp  correlation  helps  to 
illustrate  the  type  of  flow  occuring  in  the  medium. 

No  attempt  has  been  made  to  mathematically  describe 
turbulent  flow  in  porous  media  due  to  the  time  constraints 
on  this  study.  However,  this  research  has  established  the 
preliminary  work  for  further  studies  using  the  laser  doppler 
ve loci  me  ter .  It  has  also  helped  to  indicate  what  some  of  the 
desirable  features  of  future  research  should  be.  Eventually 
more  random  packs  of  different  sphere  diameter  should  be 
used.  An  attempt  was  made  to  view  the  intersecting  volume 
through  the  lucite  spheres.  Results  indicated  that 
measurements  may  be  obtained  with  the  exisiing  velocimeter 
and  a  random  pack.  There  were,  however,  refraction  and 
attenuation  of  the  laser  beams.  The  difierence  between  the 
refractive  indices  of  the  lucite  and  the  oil  was  responsible 
for  these  effects.  This  would  suggest  that  the  refractive 
indices  should  be  matched  and  a  laser  unit  of  higher  output 
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Table  7,  Comparison  of  Results  from  Cell  #1 


Reference 

( 


McFar 1  and 

This  work 


k  F  Nr p  Nr 

x  1  0"  7  m2  )  (  1  /m ) 

3.5164  70.73  0.241  123 

4.7395  48.60  0.179  104 


The  Reynolds  numbers  quoted  refer  to  the  values  at  the 
onset  of  turbulence  (  y  =  1  ). 
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Table  8.  Comparison  of  Nr  with  those  in  literature 


Author 


Pack  Description 


N 


r 


Schneebe 1 i 


27  mm  spheres  &  1  in. 
gravel  -  random 


60 


Wegner , Karabel as 
&  Hanratty 

Kyle  <S  Perrine 

Kingston  &  Nunge 
McFar 1  and 

This  work 


3  in.  sphere  -  cubic 


2-D  arrays  of  1.5  in. 
cy 1 i nders 


1.5  in.  spheres 

0.5  &  0.25  in. 
cubic  packs 


0.5  in.  cubic  pack  #1 

pack  § 2 
pack  # 2A 


90  -  120 


55  -  137 


80 


123  -  251 


104 
41  .7 
42.2 


8.  Conclusions 


This  study  has  successfully  demonstrated  the  use  of  the 
laser  doppler  velocimeter  for  flow  measurements  in  a  porous 
medium  over  a  range  of  flowrates.  As  a  result,  the  following 
conclusions  can  be  made. 

1.  A  porous  medium  should  not  be  represented  by  a  series 
of  straight  flow  channels  of  converging  and  diverging  cross 
sections  as  flow  through  such  channels  is  mainly  rectilinear 
and  approximates  that  through  straight  pipes.  The  flow 
geometry  should  include  the  curvilinear  path  of  the  fluid. 

2.  The  region  between  the  Darcy  and  the  turbulent  regimes 
may  be  described  as  a  zone  of  increasing  fluctuation  in  the 
flow.  The  extent  of  this  intermediate  zone  depends  on  the 
flow  geometry . .  F low  which  is  mainly  curvilinear  gives  rise 
to  an  intermediate  zone  much  shorter  than  that  arising  from 
rect i 1 i near  flow. 

3.  The  Forchheimer  equation  may  have  to  be  modified  if  it 
is  to  be  used  to  predict  the  pressure  drop  in  the  turbulent 

region . 

4.  The  choice  of  the  length  parameter  in  the  definition  of 
the  Reynolds  number  is  yet  to  be  established.  However  the 
Reynolds  number  as  defined  by  N  rp  suggests  that  KF^  is  a 
reasonable  choice  as  Nrp  has  been  shown  to  have  a 
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mechanistic  base. 

5.  The  laser  doppler  velocimeter  is  an  excellent 
instrument  for  non- i ntrus i ve  flow  studies.  It  can  be  used 
for  quantitative  and  qualitative  studies  of  the 
characteristics  of  flow  regimes  and  in  particular  the 
transition  to  turbulence. 

6.  With  effective  flow  diffusers,  the  location  of  the 
pressure  taps  in  the  cell  is  not  critical. 


9 .  Reccmmenda  t i ons 

The  effort  to  describe  turbulent  flow  in  a  porous  medium  is 
a  continuous  process.  To  this  end,  the  following 
recommendations  are  made. 

1.  With  the  present  velocimeter,  further  flow  studies 
should  be  performed  using  cells  with  packing  arrangements 
which  would  eliminate  rectilinear  flow.  This  would  serve  to 
indicate  how  the  onset  of  turbulence  is  affected  by  the 
curvature  of  the  flow. 

2.  Further  experiments  with  cell  2  should  be  conducted  by 
increasing  the  flowrate  well  beyond  those  encountered  in 
this  study.  This  is  to  verify  that  the  onset  of  turbulence 
indicated  in  this  study  is  indeed  so.  The  cell  has  not  been 
tested  to  see  if  it  could  withstand  the  higher  pressure 
associated  with  higher  flowrates.  This  is  an  important 
factor  when  considering  future  designs  of  cells  which  are  to 
accommodate  turbulent  flows. 

3.  A  random  packing  of  smaller  diameter  spheres  should 
eventually  be  use  to  simulate  the  porous  medium.  This  would 
require  the  use  of  a  more  powerful  laser  unit.  The 
refract i ve  i nd ices  of  the  lucite  spheres  and  the  oil  shou 1 d 
be  matched.  The  laser  beams  will  then  travel  through  the  oil 
and  the  spheres  with  mimimal  refraction.  Ati.enuation  in  the 
light  intensity  would  give  rise  to  a  lower  s i gna 1  - to-noi se 
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ratio.  It  may  then  be  necessary  to  switch  to  a  form  of 
signal  analysis  other  than  frequency  trackering  as  used  in 
this  study. 

4.  As  there  is  a  limit  to  the  gravity-fed  pressure  head 
available,  it  may  be  worthwhile  to  investigate  if  an  input 
pump  could  be  used  without  introducing  perturbation  to  the 
flow.  In  this  case,  care  must  be  taken  to  isolate  the  laser 
velocimetry  equipment  from  the  machine  vibration. 

5.  The  needle  valve  for  controlling  the  flow  should  be 
placed  as  far  upstream  as  possible  from  the  inlet  to  the 
cell  to  reduce  flow  perturbation  which  results  in 
fluctuating  pressure  changes. 
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m 

N 

f 

N 

r 

N 

rp 

Q 
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Nomenclature 

=  cell  cross-sectional  area,  cm 2 
=  constant  in  Forchheimer  equation 
=  constant  in  Forchheimer  equation 
=  friction  factor,  dimensionless 
=  constant  in  Equation  3 
=  average  particle  or  sphere  diameter,  m 
=  beam  diameter  at  1/e2  intensity,  m 
=  diameter  of  intersecting  volume,  m 
=  unit  vector  of  incident  beam 
=  unit  vector  of  scattered  beam 
=  viscous  resistance  coefficient,  m-2 
=  inertial  resistance  coefficient,  m  1 
=  focal  length  of  lens,  m 
=  doppler  frequency,  Hz 
=  incident  beam  frequency,  Hz 
=  frequency  of  measured  signal,  Hz 
=  scattered  beam  frequency,  Hz 
=  permeabi 1 i ty ,  m2 
=  fluid  mass  flowrate,  g/s 
=  constant  in  Equation  3 
=  number  of  fringes  in  measuring  volume 
=  Reynolds  No.  based  on  sphere  diameter 
=  Reynolds  No.  defined  for  porous  medium 
=  volumetric  flowrate,  cm" / s 
=  superficial  velocity,  cm/s 


. 
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local  velocity,  m/s 
time  averaged  local  velocity,  m/s 
fluctuation  of  local  velocity,  m/s 
width  of  intersecting  volume,  m 
horizontal  distance  along  cell,  m 
pressure  drop  across  taps,  Pa/m 
fringe  spacing,  m 

transit  time  of  particle  crossing  fringes,  s 

distance  between  pressure  taps,  m 

refractive  index  of  medium 

refractive  index  of  liquid 

beam  intersecting  angle,  deg 

beam  intersecting  angle  in  liquid,  deg 

laser  wavelength,  m 

laser  wavelength  in  liquid,  m 

fluid  viscosity,  Pa.s 

3 

fluid  density,  Kg/m 


. 
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Append i x  1 

Calculation  of  Properties  of  Optical  Intersection 

a)  Using  the  geometry  formed  by  the  two  intersecting 
laser  beams,  the  angle  9/2  is  calculated  from  the 
equation  given  below. 


0/2 

.  -1 
=  tan 

(  a/2 
\ b 

a 

=  3.85 

cm 

b 

=  10.20 

cm 

0/2 

.  -1 
=  tan 

l2  * 

o 

• 

o 

II 

82 


b)  Intersecting  Volume 
Given 


6/2  =  10.7° 

X  =  6328  x  10"10 

=  0.65  x  10‘3  m 

e 

f  =  0.13  m 


i  ) 


Using  Equation  22, 


w 


=  A  .  U)(f) 


IT 


D 


1 

sin  (e/2) 


4  _  (6328  x  1 0~ 1 0 )  (0, 1 3) 

71  (0.65  x  10"3)  (sin  10.7) 

-3 


=  0.869  x  10 


in 


ii)  Using  Equation  23, 

H  _  4  (A)(f)  1 

d  '  7  '  fcos  (T/  2T 

4  .  (6328  x  10~1Q) (0. 1 3) 

77  (0.65  x  10~3)  (cos  10.7) 

=  0.164  x  10"3  m 


iii)  Using  Equation  13,  the  fringe  spacing  is  given  by 


_  x 

2  si n ( e/2 ) 

6328  x  10"10 

2  sin  10.7 

=  1.70  x  10"6  m 


iv)  Using  Equation  24,  the  number  of  fringes  in  the 
intersecting  volume  is  given  by  : 


‘  fr  As 

1,64  x  10~4 

1.70  x  10'6 

ss  96 

c)  Using  Equation  16,  the  local  velocity  measured  by  the 
velocimeter  is  given  by  : 

x 

2  sin (0/2) 

(6328  x  10~10) 

2  sin  10.7 

1.70  x  10'3  ms  /kHz 


where  f^  is  measured  in  KHz. 


. 


d) 


Size  of  measuring  volume  in  relation  to  cross  section 
of  pore  space 
Consider  Cell  1 
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Area  of  widest  cross  section 
of  pore  space 

Area  of  measuring  volume 

Fraction  of  pore  occupied 
by  measuring  volume 


=  1 61 .29  mm^ 
=  1.158  mm^ 


=  1.158 
161 .29 

=  0.0347 
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Appendix  2 

Sample  Calculation  of  Flow  Data 
Consider  first  data  point  of  Cell  1 

a)  Flowrate, 


Q 


M/p 

9.0  x  10" 

850.1 

10.6  cc/s 


m3/s 


b)  Superficial  Velocity, 


q  =  Q/A 


10.6 

16.42 

=  0.645  cm/s 


c)  Reynolds  Number, 


kFbpq 

y 

(4. 7395  x  10~7) (48. 60) (850. 1 ) (0.645  x  10~2) 

16.4  x  10"3 


0.007 
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d)  Friction  factor, 


1  Ap 

—  .  •  ■  -m  ■ 


Ww^0-34  x  1q3) 

(850.1 ) ( 0 . 645  x  10'2)2 


=  198 
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Appendix  3 

Table  A.1  Experimental  Flow  Data  for  Cell  #  1 


TRANSDUCER  OUTPUT 
Volts 

FLOW 

READING 

FREQUENCY 

f  .  ,  KHz 
d 

0.11 

5 . 4  1 

21  .8 

0.23 

10.1 

65.0 

0.42 

16.  1 

139. 

0.57 

21  .0 

202. 

0.73 

21 . 0/20 . 02 

268. 

0.86 

22.8/34.0 

325. 

1  .02 

22.0/68.0 

400. 

1.15 

22.7/85.0 

440. 

1  .31 

32. 003 

495. 

1  .40 

30.30 

530. 

1  .49 

28.00 

565. 

1  .57 

27.60 

585. 

1  .65 

26.30 

610. 

1  .74 

25.10 

640. 

1  .85 

23.60 

665. 

1  .92 

23.  10 

675. 

1  .95 

22.85 

695. 

2.02 

22.35 

710. 

2.10 

21.90 

720. 

2.20 

21  .25 

730. 

2.50 

19.90 

735. 

2.80 

19.00 

760. 

3.10 

18.05 

3.31 

17.00 

3.56 

16.55 

3.94 

15.50 

4.32 

14.75 

4.65 

14.25 

Rotameter  #1 

Rotameter  #1/#2 

connected  in  parallel 

Time  in  seconds 

for  5670  cc.  to  pass 

through 

Trident  flowmeter 


. 
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Table  A. 2  Processed  Data  for  Cell  #  1 


PRESSURE  DROP 
aP/ax,  KPa/m 

Q 

cc/s 

q 

cm/s 

% 

cf 

0.34 

10.60 

0.645 

0.007 

198. 

0.68 

25.29 

1  .544 

0.018 

72.  1 

1.25 

50.35 

3.066 

0.035 

33.7 

1  .76 

73.52 

4.477 

0.051 

22.2 

2.27 

92.46 

5.631 

0.064 

18.  1 

2.61 

113.5 

6.913 

0.079 

13.8 

3.  15 

136.8 

8.332 

0.095 

11.5 

3.55 

152.3 

9.273 

0.105 

10.4 

4.06 

177.0 

10.78 

0.  123 

8.80 

4.36 

186.9 

11.38 

0.  129 

8.52 

4.67 

202.3 

12.32 

0.140 

7.70 

4.97 

205.2 

12.50 

0.142 

8.05 

5.  18 

215.4 

13.  12 

0.  149 

7.61 

5.55 

225.6 

13.74 

0.156 

7.43 

5.85 

240.0 

14.62 

0.  166 

6.93 

6.09 

245.2 

15.  14 

0.170 

6.91 

6.  19 

247.9 

15.20 

0.172 

6.87 

6.43 

253.4 

15.43 

0.  175 

6.83 

6.70 

258.6 

15.75 

0.179 

6.83 

7.00 

266.5 

16.23 

0.184 

6.72 

8.02 

284.6 

17.33 

0.  197 

6.75 

9.03 

298.  1 

18.16 

0.206 

6.93 

10.01 

313.8 

19.11 

0.217 

6.93 

10.65 

333.2 

20.29 

0.231 

6.55 

11.50 

342.2 

20.84 

0.237 

6.70 

12.82 

365.4 

22.25 

0.253 

6.55 

14.04 

384.0 

23.39 

0.266 

6.49 

15.22 

398.9 

24.29 

0.276 

6.52 
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Table  A. 3  Experimental  Flow  Data  for  Cell  #  2 


TRANSDUCER  OUTPUT 
Volts 

FLOW 

READING 

FREQUENCY 
fd,  KHz 

0.21 

2 . 3  1 

8.  1 

0.54 

5.2 

21  .8 

0.74 

6.4 

28.0 

0.92 

7.4 

33.0 

1  .  14 

8.4 

37.8 

1  .46 

9.7 

43.0 

1  .77 

11.0 

47.5 

2.07 

12.0 

48.5 

2.40 

13.1 

50.5 

2.67 

13.9 

51.0 

3.06 

15.0 

51  .5 

3.49 

15.9 

52.0 

3.75 

16.8 

51  .5 

4.04 

17.5 

51  .0 

4.44 

18.4 

50.5 

4.86 

19.3 

49.0 

5.08 

19.8 

47.0 

5.44 

20.6 

46.0 

5.70 

15. 9/24. 82 

45.5 

5.88 

21.4 

45.0 

6.04 

21 .8 

44.5 

6.50 

16.0/33.8 

7.08 

18.2/28.0 

7.54 

17.4/36.5 

8.24 

18.9/35.0 

8.95 

19.2/42.5 

1 

2 


Rotameter  #1 

Rotameter  #1/#2  connected  in  parallel 
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Table  A. 4  Processed  Data  for  Cell  #  2 


PRESSURE  DROP  Q  q  N  Cf 

A  P/ax,  KPa/m  cc/s  cm/s 


0.109 

3.29 

0.277 

0.027 

43.  1 

0.301 

9.99 

0.841 

0.082 

12.8 

0.421 

13.5 

1.14 

0.111 

9.81 

0.520 

16.5 

1  .39 

0.135 

8.16 

0.646 

21  .2 

1  .78 

0.174 

6.14 

0.862 

24.  1 

2.03 

0.198 

6.05 

1  .01 

28.5 

2.40 

0.234 

5.32 

1.18 

32.3 

2.72 

0.266 

4.79 

1  .37 

36.6 

3.08 

0.300 

4.35 

1  .53 

40.0 

3.37 

0.328 

4.08 

1  .77 

44.8 

3.77 

0.368 

3.76 

2.02 

49.4 

4.  16 

0.406 

3.53 

2.  19 

53.5 

4.51 

0.440 

3.25 

2.37 

56.9 

4.79 

0.468 

3.12 

2.61 

61.2 

5.  15 

0.502 

2.98 

2.87 

65.3 

5.50 

0.536 

2.87 

3.00 

67.3 

5.66 

0.553 

2.82 

3.22 

71.5 

6.02 

0.587 

2.68 

3.39 

73.  1 

6.  15 

0.600 

2.70 

3.50 

75.3 

6.34 

0.618 

2.63 

3.61 

77.0 

6.49 

0.633 

2.59 

3.90 

81  .2 

6.84 

0.667 

2.52 

4.28 

86.5 

7.28 

0.710 

2.44 

4.57 

90.2 

7.59 

0.741 

2.39 

5.02 

5.48 

96.0 

103.0 

8.08 

8.67 

0.788 

0.846 

2.32 

2.20 

. 
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Table  A. 5  Experimental  Flow  Data  for  Cell  #  2A 


TRANSDUCER  OUTPUT 
Volts 

FLOW 

READING 

FREQUENCY 
fd  ,  KHz 

0.24 

2 . 7  1 

10.0 

0.34 

3.8 

14.0 

0.60 

5.7 

23.0 

0.84 

7.  1 

30.0 

1  .09 

8.4 

36.5 

1  .42 

9.8 

43.0 

1.71 

10.9 

47.5 

2.  16 

12.5 

52.5 

2.92 

14.9 

57.0 

3.35 

16.0 

58.5 

3.85 

17.3 

55.0 

4.32 

18.4 

54.0 

4.90 

19.7 

53.5 

5.20 

20.4 

53.0 

5.83 

21  .6 

52.5 

6.00 

18. 0/20. 02 

52.0 

6.29 

22.5 

6.82 

16.3/37.0 

7.61 

16.2/46.5 

8.40 

15.3/58.5 

8.95 

15.6/64.0 

1 

2 


Rotameter  # 1 

Rotameter  #1/#2  connected  in  parallel 
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Table  A. 6  Processed 

Data  for 

Cel  1  U  2A 

PRESSURE  DROP 

Q 

q 

Nrn 

C  f 

A  P/A  x ,  KPa/m 

cc/s 

cm/  s 

r  p 

T 

0,13 

4.  12 

0.347 

0.034 

33.44 

0.  19 

6.47 

0.545 

0.054 

19.22 

0.34 

11.41 

0.961 

0.095 

11.28 

0.47 

15.53 

1  .31 

0.130 

8.45 

0.62 

19.53 

1  .64 

0.  163 

7.08 

0.80 

24.23 

2.04 

0.202 

5.94 

0.98 

28.23 

2.38 

0.236 

5.36 

1  .23 

34.11 

2.87 

0.285 

4.58 

1  .68 

44.58 

3.75 

0.372 

3.67 

1  .94 

49.76 

4.  19 

0.415 

3.40 

2.25 

55.88 

4.70 

0.466 

3.13 

2.54 

61.05 

5.  14 

0.510 

2.96 

2.89 

67.05 

5.64 

0.560 

2.78 

3.07 

70.58 

5 . 94 

0.589 

2.67 

3.47 

76.23 

6.42 

0.637 

2.59 

3.58 

78.50 

6.57 

0.652 

2.55 

3.77 

80.58 

6.78 

0.673 

2.52 

4.  10 

85.50 

7.20 

0.714 

2.43 

4.62 

93.00 

7.83 

0.777 

2.32 

5.32 

98.22 

8.27 

0.820 

2.30 

5.49 

104.0 

8.75 

0.868 

2.20 

■ 


MASS  FLOWRATE,  GM/S 
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FIG.  fl-1 


CflL I BRRT I  ON  CURVE  FOR  ROTAMETER  #1 


MASS  FLOWRATE,  GM/S 
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FIG 


R , 2  CRLIBRRTION  CURVE  FOR  ROTRMETER  *2 


~ 


PRESSURE  DIFF  .  .  KPfl 
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VOLTS 


FIG.  fl.3  CALIBRATION  GRAPH  OF  PRESSURE  GAUGE  #1 
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PRESSURE  DIFF  .  .  KPfl 
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VOLTS 


FID.  fl.4  CRL I BRRT I  ON  GRRPH  OF  PRESSURE  GRUGE  *2 


. 


